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Introduction 


In an atternpt to study the interaction between combustion and turbulence in supersonic 
flows, an assumed PDF has been employed^. This makes it possible to calculate the time 
average of the chemical source terms that appear in the species conservation equations as 

w s = J P(T, Yj, ..., Y u ) w s dTdY! ... dY n (1) 


where w s is the production rate of species s, P is the probability density function, T is the 
temperature and is the mass faction of species k. The probability density function is 

assumed to have the form 


PCr,Y lf ... Y n ) 


Fj(T) F 2 (Y lf Y k ) 
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where Fj is a Gaussian distribution and F 2 is the multivariate PDF developed by Girimaji 3 , r> 
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i.e., 


F 2 (Y,...Y„) = 5(1 -sin Y k k ' ! 

k=l 

where 

n 

s=lY k , c = 11 nPni/TdPk) 
k=l 
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1 


5 is the Dirac delta function, and n is the number of species present. The quantities 
appearing in the PDF can be determined in terms of the mean mass fractions and the sum of their 
variances, i.e., 

<5) 

where 

S = Z(Y k ) 2 , Q = I Y t 2 <«> 

In order to determine the averages indicated in equation (1), two transport equations, one 
for the temperature (enthalpy) variance and one for Q, are required. Model equations are 


formulated for such quantities and they can be written as 



In the above equation D is the diffusion coefficient, p is the viscosity, C t and Cq are model 
constants and t is the turbulent time scale. 

As is seen from equations (7) and (8), the turbulent time scale controls the evolution of 
both (T")^ and Q. An algebraic model similar to that used by Eklund et al. 4 was used in 
Reference 2 in an attempt to predict the recent measurements of Cheng et al. 5 Predictions were 
satisfactory before ignition but were less satisfactory after ignition. One of the reasons for this 
behavior is the inadequacy of the algebraic turbulence model employed. Because of this, the 
objective of this work is to develop a k-co model to remedy the situation. 
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Results and Discussion 

A full kinetic model based on the work of Jachimowski 6 is employed. It consists of 
eighteen reactions (Table 1 ) and allows for the presence of H? , O 2 , H, O, OH, H 2 O, HO 2 , H 2 O 2 
and treats N 2 as an inert A schematic of the experimental setup is shown in Figure 1, and 
Table 2 gives a summary of the burner exit conditions. Measurements were taken at axial 
locations of 0.85, 10.8, 21.5, 32.3, 43.1, 64.7, and 86.1 inner diameters downstream of the 
burner. 

It became evident that traditional model constants in the k-G) model were not appropriate 
for modeling the experiment under consideration. Because of this, we decided to develop, as a 
first step, a one-equation turbulence model based on the energy equation and then use it to 
develop the k-to model. This first step is now complete. 

Figures 2-4 show comparisons of mean temperature, mean H 2 O, and rms H 2 O 
distributions. As is seen from the figures, there is a need to implement the full k -10 model. 

The paper will implement the complete turbulence model and will present detailed 
comparisons with the measurements of Ref. 5. 
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Tables 

Table 1. Chemistry Model 


Reaction i 
Number 

Reaction 

A 

h 

la i; 

i 

1 

H 2 + 0 -> — OH + OH 

. 170E 4- 14 

0.00 

24157.0 

2 

H + O2 — OH + O 

. 120E 4- 18 

-0.91 

8310.5 

3 

OH + H 2 — H 2 0 4- H 

220E 4- 14 

0.00 

2591.8 

4 

C ) + //•>■ — OH + H 

506E 4- 05 

2.67 

3165.6 i 

5 

OH +OH — H 2 0 + 0 

630E 4- 13 

0.00 

548.6 j 

6 

HO 9 + H — H 2 "b O2 

.130E + 14 

0.00 

0.0 ! 

7 

HO-> + H — OH + OH 

. 1 5 0 E ~b 15 

0.00 

503.3 ! 

8 

H 0 ^ -b O — * O2 + OH 

.200E 4- 14 

0.00 

0.0 j 

9 

HO2 + 0H - h^o + o^ 

.200E 4- 14 

0.00 

0.0 ; 

10 

HO 2 + H 2 ~ HnO, 4- H 

.301 E 4- 12 

0.00 

9411.2 j 

11 

HO2 + H 0 2 — H2O2 + O2 

200E 4- 13 

0.00 

0.0 

12 

H + H 2 0 2 — H2O + OH 

.100E 4- 14 

0.00 

1801.7 

13 

0+ H2O2 — HCh + OH 

.280E 4- 14 

0.00 

3220.9 

! 14 

OH + H2O2 — H ?0 4- HO2 

.700E 4- 13 

0.00 

722.2 

1 15 

H +OH + \f — H 2 0 4- M 

.221E 4- 23 

-2.00 

0.0 

! 16 

H 4- H + M — * H 2 4- A/ 

.730E + 18 

-1.00 

0.0 : 

1 17 

H + O2 + M — HO2 + M 

230E + 19 

-1.00 

0.0 ! 

18 

H2O2 -r M — OH + OH 4- M 

121 E 4- 18 

0.00 

22898.8 ij 


Jnits of *4 are a multiple of cm 3 mole' 


Table 2. Exit Conditions 


lauie 

Exit Conditions 

Hydrogen Jet | 

Outer Jet 

Ambient Air j 

Mach Number 

1.0 

2.0 

r— 0. —\ 

Temperature. K 

545 

1250 

CO 

O 

O 

Velocity, m/s 

1780 i 

1417 

0. 

Pressure. MPa 

■112 

.107 

101 

Mass Fraction 


Yfii 

1. 

0. 

0. 

Yo , 

0. 

.245 

233 

— *1 

Ysi 

0. 

.58 

757 



ii — 1 — 

Yh-,0 

0. 

.175 

ij.01 
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Figure 1: Schematic of test apparatus 
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Figure 2: Mean Temperature Comparison With Experiment 









